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Abstract. We establish bounds on the energy of a system of N identical 
bosons bound by attractive pair potentials and obeying the semirelativistic 
Salpeter equation. The lower bound is provided by a 'reduction', with the 
aid of Jacobi relative coordinates, to a suitably scaled one-body Klein-Gordon 
problem. Complementary upper energy bounds are provided by means of a 
Gaussian trial function. Detailed results are presented for the exponential pair 
potential V{r) = —vexp[—r/a). 
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1. Introduction 

We consider a system of N identical bosons which interact by attractive central pair 
potentials V{rij). There is no restriction in general concerning the number d of spatial 
dimensions, but, for definiteness in the presentation, we shall assume d — Z. We 
suppose that the potential has the form V{r) = vf[r/a), where r = |r|, w > is the 
coupling parameter, f{r/a) < is the shape of the potential, and a > is a range 
parameter. The Salpeter semirelativistic Hamiltonian for a system of N identical 
bosons of mass m and individual momenta pj is given by 

N N 

i—l l—i<^j 

where r^j = jr.; — rj\ and pi — \p^\. Our principal goal is to find or estimate E, 
the lowest energy of this system. The ease with which the A^-body problem can be 
posed in the Salpeter theory is the main reason for our adopting the semirelativistic 
approach. For suitable potentials, the Hamiltonian is bounded below and to a certain 
extent the spectral theory mimics the corresponding Schrodinger theory [1-3]. 

One path to a lower bound is to exploit the boson permutation symmetry of the 
exact normalized ground state ^' to effect a 'reduction' generated by the equation 
E — (^f, H'i') — ('5, 7i^), where H is the Hamiltonian for a specially scaled one-body 
problem. A lower bound to E is then given by the bottom of the spectrum of Ti.. 
Such reductions have a long history which is reviewed briefly in Ref. [4]. Different 
reductions are possible and lead to a variety of lower bounds. For example, we obtain 
immediately from Eq.(IT]) 

= = (1',i?2*), (2) 

where 



A^ 



pI + m2 + Jp^ + m2 + (A^ - l)V{ri2) 



(3) 



The reduction is not yet complete since H2 is the Hamiltonian for a scaled two-body 
problem. If we adopt new coordinates R = ri + r2 and r — ri — r2 for this problem, 
wc find Pi = P + p and P2 — P — P- Thus if we write \/pf + rn? — \p^ -\-mk\, i = 1,2, 
where fe is a unit vector orthogonal to Pi and P2, then an application of the triangle 
inequality yields 

2\p + mk\ = \(p — P + mk) + {p + P + mk)\ < \ — P2 + "^'^1 + \Pi + mk\. 
Thus we have the spectral inequalities 



H>H2> nN/2 = N^p-' + m^ + jvfir/a) > £n/2, (4) 

where 7 = ^N{N— 1) and £n/2 is the bottom of the spectrum of Hn/2 and provides a 
lower bound to E, the ground-state energy of H. We note that this use of the triangle 
inequality is equivalent to setting the centre-of-mass momentum P of the two-body 
Hamiltonian H2 to zero. 

For boson systems there may be an energy advantage (leading to a higher 
lower bound) if, from the outset, Jacobi (orthogonal) relative coordinates can be 
accommodated by the reduction scheme. Intuitively, the use of relative coordinates 
allows us to remove the centre-of-mass kinetic energy of the whole system, rather 
than dropping the centre-of-mass energy only at the final stage, after the reduction to 
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a scaled two-body problem such as One way of achieving this for the Salpeter 
problem involves squaring the reduced one-body Salpeter equation of the form 
^/jp^+m^tp = [E — V)tp. This allows us to invoke some remarkable expectation 
identities for Jacobi coordinates: these identities, which we shall refer to explicitly 
below, are proved in Appendix A. 

Since we use relative coordinates both for the upper and lower bounds, in section 2 
we define the Jacobi coordinates {Pj}^2 ^^'^ ^'^^^ some key properties. An upper 
energy bound Eu = ($, H$) is provided by an application of a normalized Gaussian 
trial function This latter calculation is greatly simplified by the boson permutation 
symmetry of $ and also by the fact that this particular trial function is not only 
symmetric under the permutation of the relative coordinates but also factors in the 
form $ = Y\f^2^iPi)- upper bound is discussed in section 3. In section 4 we 

introduce a variant of the Jacobi system suitable for the derivation of a lower bound 
expressed in terms of a Klein-Gordon one-particle problem. The proof assumes that 
this Klein-Gordon problem has a discrete eigenvalue that is either (i) positive or (ii) 
monotone increasing as a function of the mass. In section 5 we consider the example 
of the attractive exponential potential. 

2. Relative coordinates 

Since boson permutation symmetry is naturally expressed in terms of the individual- 
particle coordinates {r^l^^, it might seem an unwelcome complication to introduce 
relative coordinates {Pi}fL2- Their advantage arises when we consider the reduction 
(H) = 1(H), where H is a scaled one-body problem. Jacobi relative coordinates may 
be defined with the aid of an orthogonal matrix B transforming old {r^} to new 
{pj} coordinates by [p] ~ B\r]. The first row of B, with all entries Bu = 
defines a centre-of-mass variable p^, its second a pair distance P2 = (t"i — r2)/V^, 
and the fcth row (fc > 2) first has k — 1 entries Bu = l/\/k{k — 1), the fcth entry 
Bkk = —a/ {k — l)/k, and zero for all remaining entries. The momenta {TTj} conjugate 
to the {pj} read [tt] = {B~^)'^[p] = B[p]. Now, for an A''-boson problem with an 
attractive potential V{r), let ^(p2, p^, - ■ ■ , p^) be the exact normalized ground-state 
eigenfunction corresponding to the lowest energy E. Boson permutation symmetry is 
a powerful constraint that greatly reduces the complexity of this problem. Although 
a non-Gaussian wave function is not necessarily symmetric in the Jacobi coordinates, 
we do have the remarkable AT-representability expectation expressions from Appendix 
A: for i,j>l 

(*,(p,.p,.)*) =5,,(*,p2vl/), (5) 
{^,{7Vi-7Tj)^) =%(*,7r2*). (6) 

3. A Gaussian upper bound 

We now take advantage of the special properties of the Gaussian trial function $ to 
effect an upper reduction of (H) . In the centre-of-mass frame of the A''-body system, 
where tti = 0, we have from the last row of the matrix B 

Pj^ = -L=7ri - VXttn = -a/Attjv, A = ^ ^ . (7) 
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Meanwhile our trial wave function $ is not only a boson function but it is also 
symmetric and factors in the relative coordinates. We have explicitly: 

N , , I 

where a > is a variational scale parameter. The boson permutation symmetry of 
the trial function allows us to write E < Eg = ($, , where we have 



1^$, N^\T,l + m^+^V{\V2p^\) $^ 



Setting r = ri — r2 = \f^P2 and P = ^^2/^, we find E < Eg — {(j),Hg(j)), where 

Hg = 7VV2Ap2 + m2 + jV{r). (8) 

We now introduce some new parameters which will allow us to express both the 
upper and lower energy bounds in compatible succinct forms. With V{r) ~ vf{r/a), 
X={N -l)/N and 7 = \N{N - 1), we now define 

^ = ma/(2A)* e^Ea/{2-f^), Cg ^ Ego/ i' = "f^va/2. 
The Gaussian upper bound then becomes 



e < Cg — Tii\Q.{(l),hg(j)) ^ where hg = \/pP^+Jj? + v f{r). (9) 

(T>0 

We note parenthetically here that, if the same scaling is used for the reduced 
Hamiltonian 7ijv/2 of El- O? the result h2 is comparable with hg and is given by 



h2^^W/W)+7+^f{r), \^{N~l)/N. (10) 

We shall use this later in the paper to show that the N/2 bound is weak: we shall 
consider 

e2g = min(</), /i20) > 62, (11) 

o->0 

where 62 is the bottom of the spectrum of h2, and observe (in particular cases) that 
&2g < Cfe, where ek is the Klein-Gordon lower bound that we shall derive in the next 
section; hence, for these problems, we know, without calculating 62, that 62 < e^. 

However, our main point here is that the upper energy bound can be written 
in terms of an optimized expectation (with respect to a Gaussian wave function) 
of a special one-particle Salpeter problem with Hamiltonian hg. The kinetic-energy 
expectation value may be expressed in terms of the modified Bessel function of the 
second kind Ki{x) [7]. Evaluating the integral, for convenience, in momentum space, 
yields 



=^/(x), 
V 2x 



where a = and the integral I{x) is defined by 

I{x) = — y/2x + t^e-'*\^dt = —xe^'Kiix). 
Vi" Jo V'^ 

Meanwhile, the potential-energy expectation becomes 

(0,/(r)^) = j'^^ 
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where the integral J{y) is given by 



J{y) = -r 



-f{yt)t^dt. 



We now use s — l/a = as a variational parameter and we obtain 

,2\ 



e < e„ 



s>0 



si 



2s2 



(12) 



4. The lower energy bound 

In order to effect a new reduction that incorporates tti = 0, we first choose to use ppf 
and Ptv-i ^'^d obtain the reduction 



N 



N 



+ YVPjv-i +7V'(|rAr_i - rAr|) 



where 7 = iA^(A^— 1). We note that the last column of B (which defines tt^v) is given 

by 

[b,b,b,...,b,-{N -l)b], where 6= l/^/N{N-l). 

Now we express the equation for E in terms of new coordinates defined by the following 
relations 



R 

r 



P a 
a —(3 



Pn 

Pn-1 



P 

P 



where 



N -1 



Consequently, 



and 



^ = A 



f3 



{N-l)b = a~^, Nb^a, S ^ I + 5^ = 2X. 



r = |r|, 




(3 a 




a -f3 




\l N 


2 






N-1 




N ' 




1 + 5^ =2X 



< 1. 



We note that because of the boson permutation symmetry, and after we set the centre- 
of-mass momentum variable tti = 0, we have generally that 

{<\Pn-i\)) - {<\P-5P\)) = (ni\p + SP\)) = («.(|p^|)), 

where K{p) is any appropriate kinetic-energy function. The expression for the energy 
now has the form E = {Ti) , where 



n = N^{p + spy + m2 + -fV{r). 



(13) 



The Hamiltonian Ti. is bounded below by the simple bound TiN/i of Eq. This 
result is a consequence of the inequality 



2(|p + mfc|) ^ {\{p ~ 5P + mk) + {p + 5P + mk)\) < 2{\p + dP + mk\ 
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We now consider the eigenequation Tii/i = £■(/;, where ■ip{r,R) E V C L^{R^) retains 
some of the boson permutation symmetry imphcations of the full wave function, 
namely from Eq. (p^) = (-P^), (p ■ P) = 0, and hence 

{{p + 6Pf) = {p^ + d'P') = {l + S^){p') = 2A(p2). (14) 

Let £ be the minimum of {ip, Ti-tlj) corresponding to normalized ip satisfying (|14|) and 
let be the exact normalized iV-body wave function, then 

(*,i7*) = (5',^^') > i^,ni;) = £. 

The eigenvalue equation for £ reads 



N^/{p + spy + to2 ^ = (£ - -/vf{r/a)) ^. (15) 

Squaring these equal vectors and use of (fT4|) . and also the scaling change (to new 
dimensionless variables) {r,p} {ra,p/a}, implies 

> (/ + i^/(.) - I^/^(r)) ^) . (16) 

In contrast to our upper bound, where we worked directly with the Salpeter 
Hamiltonian Q, for the lower bound we are led to consider a complementary lower 
Klein-Gordon operator. If the minimization on the right-hand side of (fT6|) could be 
effected exactly (as it can for the gravitational problem f{r) = —1/r [5]), then it would 
remain to extract a lower bound to £ from the resulting inequality. For the general 
problem that we are considering here, we have to find another approach. To this end 
we have explored what happens if we set 5 = £fc in Eq. (fT6|l . and solve the resulting 
Klein-Gordon equation for the ground-state energy e. This Klein-Gordon equation is 
given explicitly by 

{{£kaf~N^{maf)(t> 

-4,f,^+(^^/(.)-3l!;;^/^(.)V. (17) 



This exercise is only useful if an a priori relation can be established between £k 
and £. More particularly, we shall now establish a relationship between £k and 
the lower bound to £ provided by the inequality given in Eq. (|16p . We employ a 
similar parameter set to that for the upper bound, that is to say A = {N — l)/N and 
J = lN{N - 1), and also 

^ = TOa/(2A)^ ek = £ka/{2j^), (18) 
We thus arrive at the Klein-Gordon equation in a standard form, namely 

{p'+fi^)cj,= {i.f{r)-ekf<j>, (19) 
and, equivalently, in Schrodinger form 

(p2 + 2efei//(r) - (ufir))^) (f> = F(efe)0 = {el - ^l^)c^. (20) 

The spectral function F{ek) in Eq. (I^Dl) describes, for each fixed coupling v, how the 
lowest eigenvalue of the Schrodinger operator p^ + 2ekvf(r) — (vf(r))'^ depends on 
the energy parameter e^; the Klein-Gordon eigenvalue exists if there is a solution of 
the equation F{ek) = — m^; these curves are shown in Fig. [T] for the case of an 
attractive exponential potential. Similarly, the original lower-bound inequality (fTH)) 
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kgsro(num): V(r) = -exp(-r) 




-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 I 



Figure 1. The curves of Schrodinger energies F{e) for v = 1 (0.25) 5.5, and 

— m? for m = 0.1(0.1)2.1 for the exponential pair potential V{r) = —ve~^. 
The intersection points yield the Klein— Gordon ground-state energies e = e{m, v). 
For these solutions to exist, the Schrodinger curves must be steeper than the 

— curves: thus we have F'{e) < 2e at the intersection points, even when 
e is negative; this is a sufficient condition for the Klein-Gordon energy e to be a 
lower bound to the corresponding Salpeter energy E > e. 

becomes an equality in terms of the parameters (fT5)) and also a new mass parameter 
p. > fj. used here as a 'slack variable'. Thus, we have: 

{p^ + 2evf{r) - {vf{r)f) <t> = F{e)<t> = (e^ - /I^)^. (21) 

In earlier work [6] on the semirelativistic and Klein-Gordon one-body problems, 
we have proved that, under appropriate assumptions, the Klein-Gordon energy is 
monotone in the mass. More specifically, in our application, we can state this result in 
terms of Eq. (pij) . namely, for each fixed if either e > 0, or e < and e > \F'(e), then 
e'(/[t) > 0. Thus, under these conditions (which must be verified for each potential), we 
have from Eqs. ([^0)1 and (pTjl that e > e^. Thus £fe provides us with a lower bound to 
the A^-body ground-state energy E. In terms of the dimensionless energy parameters, 
our principal results are summarized by the following inequalities: 

ekifJ^i^) < e < eg{fi,iy), (22) 
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where is given by solving the Klein-Gordon equation (j20p . Cg is given by the 
variational equation (fT^ . and the parameters are defined in equation (fT5)) . Thus, 
whenever a solution exists, we necessarily have F'(e) < 2e. For example, in the case 
of the exponential pair potential we show this relationship explicitly in Fig. [1] 

5. The exponential potential 

We now consider the example with potential shape in dimensionless form given by 
/(r) = —e~^. The lower bound ei^{fi,i>) is given by solving the explicit one-particle 
Klein-Gordon equation 

(p2+^2)^^ (g^^^g-.)2^ (23) 

and the upper bound eg(/i, ly) is given by 



where 




and 

j^y)^^JL / e-'\-yH''dt 

Jo 

= ^-^(2 + y')(l-erf(y/2))e^'/^ 

For the choice /i = 1 we exhibit some numerical results in Table 1 and Fig. [51 We 
note that the mass parameter fj, decreases 'slowly' with increasing N >2 and satisfies 
ma > /i > ma /V2. In Table 1 we include values of e2g calculated from Eq. (fTTj) for 
A'^ — > c», or A — *■ 1: since e2g < efc, we know that the Klein-Gordon bound Ck is better 
than the N/2 bound, that is to say 62 < e2g < Ck- Moreover, we see from Table 1 
that for large N , the Klein-Gordon lower bound is in fact much better than the N /2 
bound. 

The ultrarelativistic limit m — > 0, which is perhaps meaningful in the spinless 
Salpeter theory, has no useful counterpart expressed by the Klein-Gordon equation. 
The Salpeter Hamiltonian in this limit (and for the exponential example) becomes 
H = \p\ —ve~'^ : this problem is not very different qualitatively from the corresponding 
Schrodinger problem H = p^ — z/e"'', which also has a discrete eigenvalue for every 
ly > not too small. By contrast, the Klein-Gordon equation ([55)1 with /i = implies 
Ck = 0, and leads to the Schrodinger problem (p^ — i''^e~'^'^)ip = 0, which has a solution 
only for very particular v not too small. Clearly, the Klein-Gordon lower bound is 
not useful for the ultrarelativistic limit. 

6. Supercriticality and system size 

We consider first a one-body problem with an attractive (vector) potential V{r) = 
vf{r) in three spatial dimensions. We contrast four theories: Dirac, Klein- 
Gordon, Salpeter, and Schrodinger. For the singular Coulomb or gravitational case 
/(r) = — all four theories agree that binding is obtained for arbitrarily small v. 
Meanwhile, the three relativistic theories agree that the coupling v cannot be very 
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large: we have for Dirac u < 1, for Klein-Gordon v < ^, and for Salpeter v < 2/it. 
For these three relativistic theories, the discrete eigenvalues are all positive and the 
supercritical phenomenon that we discuss below does not arise. None-the-less, the 
coupling restrictions lead to limitations on the system size, as we have discussed earlier 
for a static-gravity boson-star model [4] . 

We now set aside the Coulomb case and consider bounded short-range potentials 
such as the exponential potential /(r) = — e""" discussed in the previous section. 
For this class of problems, the four theories now agree that (in three dimensions) 
binding is obtained only for v sufficiently large. The term 'critical coupling' is used 
for the coupling value Vc, in each theory, beyond which a discrete eigenvalue appears. 
As V is increased further, this eigenvalue decreases. In the fully relativistic theories 
another critical value v = Vg is reached at which E = —m. This situation is called 
'supercriticality' and Vg is the corresponding supercritical coupling. The relativistic 
Dirac and Klein-Gordon theories will not accommodate larger couplings than their 
particular supercritical values. By contrast, for the Salpeter and Schrodinger cases, 
as V is increased, the eigenvalues simply decrease without lower limit. It is now clear 
that if we use the Klein-Gordon eigenvalue to bound the Salpeter eigenvalue below, a 
limitation is introduced which was not originally present in the Salpeter problem itself, 
namely that the coupling should not be very large. One could perhaps argue that, just 
as the restriction E > —m, which emerges internally and technically from the theory 
of the Klein- Gordon equation, is 'understood' in terms of another theory, namely pair 
creation in quantum field theory; thus, significance of this constraint could be acquired 
by the Salpeter theory from the Klein-Gordon theory. All we are saying here is that, 
if the Salpeter theory claims to be partially relativistic, then it is natural to impose 



Table 1. Energy bounds < e < Cg for the A''-boson system bound by an 
attractive pair potential with shape f{r) = — exp(— r). The parameters are given 
hy ^l = ma/{2\)h = 1, e = Ea/{2'yi), u = 'yhva/2, \ = (N - 1)/N, and 
7 = iJV(JV — 1). The energy e2g is an upper bound to the N/2 lower bound 62- 
The table shows e2g for AT — ► oo. Since e2g < e^, we know that the Klein-Gordon 
bound Cfc is better than the N/2 bound, that is to say, for large N, e^. < eig < Sk- 



V 


eig 


efe 


eg 


0.6 


0.993207 






0.8 


0.962091 


0.996825 




1 


0.917724 


0.980384 


1.002270 


1.1 


0.895169 


0.967625 


0.992897 


l.r, 


().7()(il2G 


0.891979 


0.928149 


2 


0.573703 


0.754224 


0.802883 


2.5 


0.352544 


0.580890 


0.642092 


3 


0.109627 


0.380273 


0.454366 


3.5 


-0.150549 


0.157791 


0.245264 


4 


-0.424867 


-0.082840 


0.018650 


5 


-1.007410 


-0.609771 


-0.476508 


5.67 




-0.993110 


-0.833505 


5.8 






-0.904766 


6 






-1.015560 
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Figure 2. Lower (L) and upper (U) energy bounds to the ground-state energy e of 
the A'^-boson system with pair-potential shape /(r) = — exp(— r). Dimensionless 
units /i = ma/(2\)2 = 1, e = Ea/{2f2')^ v = are used, in which 

A = (N — l)/N, and 7 = ^N(N — 1). The lower and upper data are given, 
respectively, by Eqs. l(23)l and l(24)l. 

the condition that E > — m, or, what is sufficient, that v not be too large. Using the 
Klein-Gordon theory to bound Salpeter below achieves this automatically. In terms 
of the A'^-body reductions discussed in this paper, this means that the scaled coupling 
1/ — 72-ya/2 should not be too large. Since 7 — ^N{N — 1), this means, for a given 
pair interaction, that the particle number N, in turn, should not be too large. Thus 
we arrive at the conclusion that such relativistic systems in their ground states will 
not become giant balls as N is increased, but will divide into smaller aggregates. For 
the exponential example with m — 1, a — \/2, v small and N large so that A « 1, we 
have from Table 1, v = ^^vajl < 5.7, which is to say, N < llA/v. Thus we conclude 
that a phenomenon which seemed to be a special feature of the — 1/r potential may 
be more widely present in relativistic many-body systems. 
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7. Conclusion 

The principal disadvantage of the semirelativistic Salpeter equation is that it is 
not fully relativistic. It perfectly well accommodates the harmonic oscillator, a 
strong interaction not at all tolerated (as a vector potential) by the Klein-Gordon or 
Dirac theories. Moreover, the Salpeter Hamiltonian is non-local, and this introduces 
technical problems for spectral estimation. The main advantages of the Salpeter theory 
are (i) the ease with which the relativistic many-body problem can be formulated and 
(ii) the fact that the energy spectrum is characterized, and may therefore be estimated, 
variationally. For a wide class of attractive pair potentials that vanish at infinity, 
we have shown in this paper that Jacobi relative coordinates and the permutation 
symmetry of the boson states allow us to construct ground-state energy bounds in 
terms of certain reduced one-body problems: the lower bound is provided by the 
ground-state energy of a Klein-Gordon problem; and the upper bound is given by 
the expectation of a Salpeter Hamiltonian with respect to a scale-optimized Gaussian 
trial function. These bounds may well have theoretical as well as computational 
applications. For example, as N increases, the effective potential coupling v in the 
scaled one-body Klein-Gordon problem p9|) increases by the factor N{N — 1) ~ N 
whereas the kinetic-energy term remains unchanged. Even for bounded potentials, 
such as the exponential potential, this implies that this lower-bound Klein-Gordon 
problem will, for sufficiently large iV, become supercritical, with e < — /x. This may, 
in turn, indicate that the iV-body problem will not be stable for large aggregates, 
and will necessarily break up into smaller fragments. This would be the expression 
of a truly relativistic many-body phenomenon with no pro-image in non-relativistic 
theories. 
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Appendix A 

For definiteness, we consider the relative momenta tt^ and iVj with i,j > 1. Since 
[tt] = B[p\, we have 

^'•^J = ^'fePfc ] ^^^^'^ ) = XI B^kBjk (Pfc) +^ B.kBji (pfc • p;) . (Al) 

\ k / \ k ) k k=f^l 

By using the boson symmetry of the wave function we have 

\ k ) \k+l ) 
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- = ( E ( E ) " ( E B^^^i^^ ) + I E B^^^ii ■ (A3) 

k \ k / \ I / \ k / \k^l J 

li i = j, then the orthogonality of the matrix B tells us that 

E ^ ^ = - y^^BikBg. 

k k^l 

Meanwhile, if z 7^ j, we know that the corresponding rows of B are orthogonal; hence, 
in this case, J^k BikBjk — 0; thus, both of the coefficients in (A2) are zero. Hence we 
conclude for all z, j > 1 

(tt, • Tz.j) = % {nl) = ((p2) - (pi . P2)) . (A4) 

The proof for the relative coordinates {p^} is identical. 

References 

[1] E. E. Salpeter and H. A. Bethe, Phys. Rev. 84, 1232 (1951). 
[2] E. E. Salpeter, Phys. Rev. 87, 328 (1952). 

[3] E. H. Lieb and M. Loss, Analysis (American Mathematical Society, New York, 1996). 
[4] R. L. Hall and W. Lucha, J. Phys. A: Math. Gen. 41, 355202 (2008). 
[5] R. L. Hall and W. Lucha, J. Phys. A: Math. Gen. 39, 11531 (2006). 
[6] R. L. Hall and W. Lucha, Preprint, arXiv:0906.2001 [math-ph]. 

[7] Handbook of Mathematical Functions, edited by M. Abramowitz and L A. Stegun (Dover, New 
York, 1972). 



